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Abstract 

Adopting the 3 neutrino framework, we present an updated determination of the 
oscillation parameters. We perform a global analysis and develop simple arguments 
that give essentially the same result. We also discuss determinations of solar neutrino 
fluxes, capabilities of future experiments, tests of CPT, implications for neutrino-less 
double-/? decay, (3 decay, cosmology. 
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(4.6 10) meV 

(51 66 ) meV 


(12H-57)meV 
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Table 1: Summary of present information on neutrino masses and mixings from oseillation data (upper 
rows) and inferences and limits on non-oscillation probes (lower rows) A 99% G.L. range is a 2.58a 
range. Some 0z^2/3 data are controversial, and h ~ 1 parameterizes uncertain nuclear matrix elements. 
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Figure 1: Summary of present information on neutrino masses and mixings from oscillations. 
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The most plausible extension of the Standard Model that allows to interpret a wealth of neutrino 
data piiisiiiisiiniiziEiisiiniiii] consists in adding a Majorana mass term for neutrinos, 

^ = ^sM + ■ m ■ 1 / + h.c.) m = V* ■ diag(mi, m 2 m 3 . (1) 

For normal mass hierarchy we order the neutrino masses mi^ 2,3 > 0 as mi < m 2 < m 3 , whereas for 
inverted mass hierarchy we choose m 3 < mi < m 2 . The neutrino mixing matrix is 

F = i?23(023) diag(l,e'‘^,l) RisiOis) Ruie^). ( 2 ) 

A similar result holds in the case of Dirac masses, with the difference that the number of physical 
parameters decreases from 9 to 7: the Majorana phases a and (5 can be reabsorbed by field redefinitions. 

Data on neutrino oscillations fix 612 , 623 , ^"ii 2 and |Am 23 | where Am|^- = mj — mf. As discussed 
later (sections ^ and [ 21 ) our present knowledge of oscillation parameters is approximatively summarized 
in the upper rows of table Q The uncertainties are almost Gaussian in the chosen variables; hg. ^ 
shows the full functions. Correlations among parameters are ignored because negligible, with the 
exception that the upper bound on 613 depends on |Am| 3 |. While 612 is rather precisely measured, 
the other 2 mixing angles have large uncertainties. Planned long-baseline oscillation experiments can 
strongly improve on Am 23 and (if 9\3 > 1 °) measure @ 13 , the phase (p, and determine which type of 
mass hierarchy is realized in nature. 

Oscillation experiments, however, are insensitive to the absolute neutrino mass scale (say, the mass 
of the lightest neutrino) and to the 2 Majorana phases a and (3. Other types of experiments can study 
some of these quantities and the nature of neutrino masses. They are: /3-decay experiments, that in 
good approximation probe mPp^ = (m ■ m^)ee = i neutrino-less double-beta decay {^u2f3) 

experiments, that probe the absolute value of the Majorana mass mee = Si cosmological 

observations (Large Scale Structures and anisotropies in the Cosmic Microwave Background), that in 
good approximation probe rucosmo = • 93.5eV = Si^*- Only neutrinoless double beta decay 

(0i/2/3) probes the Majorana nature of the mass. The values |m,ee|,rucosmo are unknown, but 
can be partially inferred from oscillation data. Tabled shows our results, discussed in section (21 in 
the limit where the mass of the lightest neutrino is negligible. In the opposite limit neutrinos are 
quasi-degenerate and |mee|,mcosmo can be arbitrarily large. 

From the point of view of 3 massive neutrinos, it is natural to divide in three parts a discussion of 
the present situation and of the perspectives of improvements, namely: 

1. Oscillations with ‘solar’ frequency, that tell Am ^2 ^rid 612 and give a sub-dominant constraint on 
O 13 . In section n we discuss solar and reactor neutrino experiments, showing that the program 
of measurement of parameters is well under way (if not accomplished), and discussing other 
interesting related goals. 

2 . Oscillations with ‘atmospheric’ frequency, that tell Am 23 , ^23 and 613, are discussed in section|21 

3. Non-oscillation experiments, that can tell the absolute neutrino mass. In section Owe discuss 
the present status and assess the implications of the existing information on neutrino oscillations 
for these experiments, particularly for 0i/2/3. We conclude by commenting on the recent claim 
of evidence for this transition iniiTB]. 

We assume the 3 neutrino framework because it is plausible, well defined, restrictive and compatible 
with data. However it is just an assumption, and before proceeding we recall some alternatives. The 
most plausible one is the presence of extra light fermions (‘sterile neutrinos’) or bosons, which can 
manifest in many ways. Going to rather exotic scenarios, Lorentz or CPT invariance (here studied 
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in fig. n might be violated in neutrinos, that might have anomalous interactions (gauge couplings, 
or magnetic moments, or else), might not obey the Pauli principle, etc, etc, etc. Present solar and 
atmospheric data cannot be explained by these alternatives, which however might be present as sub¬ 
leading effects on top of oscillations among active neutrinos, such that our determinations of active 
oscillation parameters would need model-dependent modifications. The present data do not give 
any clear indication for extra effects, but contain some anomalous hints. Most notably, the LSND 
anomaly M is not compatible with the 3 neutrino context we assume. 


1 Oscillations with solar frequency 

A few years ago the solar anomaly rested on global fits that combined solar model predictions with a 
few measurements of solar neutrino rates. In recent times, the situation changed. As prospected in m 
(written a few years ago, while sub-MeV solar experiments were discussed as a tool for discriminating 
LMA from LOW, SMA, QVO,...), SNO, KamLAND, Borexino had in any case the capability to 
identify the true solution of the solar anomaly and make precision measurements of the oscillation 
parameters. This is where we are now. From the point of view of the determination of the oscillation 
parameters, solar neutrino experiments are in a more advanced stage than atmospheric experiments, 
as clear from hg. Q KamLAND and SNO play the key role in the determination for Am ^2 612 
respectively, and almost achieved the 2.5% and 10% accuracy prospected in m- In section o we 
present a global solar fit, and in section W?I\ we show that it is dominated by very simple and robust 
inputs. In section 01 we extract from data the survival probability of low-energy neutrinos, and use 
it to study how solar data mildly restrict @ 13 . In section Fl.dl we reassess the interest in proceeding 
with low-energy solar neutrino experiments. In section 11.51 we discuss how well present solar data 
determine solar neutrino fluxes and discuss the impact of Borexino. In all cases we extract from 
simple arguments the main general resnlts, that we compare with ‘exact’ results of global fits performed 
in specific cases. 

1.1 Updated fit of solar and reactor neutrino data 

We begin by presenting a global fit of solar and reactor neutrino data assuming oscillations among 
active neutrinos with negligible @ 13 . We include 

• The final SNO CC, NC, ES rates measured during day and during night without [H] and with [7j 
salt, that gives enhanced NC sensitivity. 

• The Super-Kamiokande ES spectra [H]. 

• The Gallium rate, iZca = (68.1 ± 3.7) SNU, obtained averaging the most recent SAGE data with 
the final Gallex and GNO data ISIHli- 

• The Chlorine rate P, i?ci = (2.56 ± 0.23) SNU. 

• The KamLAND reactor anti-nentrino data with prompt energy higher than 2.6 MeV jHj. 

Solar model predictions and uncertainties m are revised including the recent measurement of the 
^^N(p, 7)^^0 nuclear cross section which reduces the predicted CNO fluxes by roughly 50%. 

The result of our oscillation fit is shown in fig. 12^, where the best fit point is marked with a dot. 
The Icr (i.e. 68 % C.L.) and the 99% C.L. (i.e. 2.58 it) ranges for the single parameters (1 dof) are 
summarized in tabled The total evidence for an effect is now about 12cj in solar v data and about 
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Figure 2; Best-fit regions at 90, 99 and 99.73% CL. Fig. assumes CPT invariance and combines 
solar V data (dashed red contours) with reactor v data (dotted blue contours). In fig. we show how 
data determine the high- and low-energy limits of P{V(, —> ufij, as precisely described in the text, 

6 cj in KamLAND i> data. Fig. EK also shows separately the fit of solar data (dashed red contours) and 
the fit of KamLAND data (dotted blue contours). 

We note in passing that, according to the Pearson yfi goodness-of-fit (gof) employed by many old 
global analyses, today the LOW solution (with Amf 2 ~ 10“'^eV^) still has a good gof. This does not 
mean that LOW is compatible with data, but happens because the Pearson test is not a statistically 
powerful gof test when dof 3> 1 m- global fits of solar data input roughly 100 experimental inputs. 
As we now discuss, the output depend almost only on a few pieces of data. 

1.2 The meaning of ‘global fits’ 

Our results are based on a careful global fit. We now point out how a simple approximate analysis is 
sufficient to get results practically equivalent to those of global fits. 

The solar mass splitting Am ^2 Is directly determined by the position of the oscillation dips at 
KamLAND, with negligible contribution from solar experiments. (More precisely, this will be rigor¬ 
ously true in the future. For the moment solar data are needed to eliminate spurious solutions mildly 
disfavored by KamLAND data, as illustrated in fig. EK- Solar data have practically no impact around 
the global minimum of the and consequently on the determination of Amf 2 as reported in table ^). 
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atmospheric oscillations 


solar oscillations 



Figure 3: Test of CPT-violating neutrino masses. We show the separate fit for Am^ in neutrinos and 
Afh'^ in anti-neutrinos, marginalized with respect to the mixing angles 612 and 623 . The atmospheric 
fit includes data from SK, K2K, Maero. The solar fit includes data from SNO, SK, Gallex, Sage, 
Homestake, KamLAND. These plots update the original results in fig.s 5 and 6 of m 


The solar mixing angle is directly determined by SNO measurements of NC and CC solar Boron 
rates. Assuming flavour conversions among active neutrinos, SNO implies^ 

(P(z/e ^ Z^e)) = = 0.357 ± 0.030. 

This should be compared with the theoretical prediction for {P{TZe ^e))-, given by a simple ex¬ 
pressions that does not depend on the solar density profile because LMA oscillations are almost 
completely adiabatic. At MeV matter effects dominate such that t'e produced around the 

center of the sun coincide with the z ^2 eigenstate in matter and exit as the V 2 eigenstate in vacuum, 
so that P{h'e —> Ve) — siu^ 012. In the energy range explored by SNO, matter effects at the production 
region are not fully dominant, such that the above approximation gets slightly corrected to^ 

(^(t'e —^ i^e)) ~ 1.15 sin^ 012 so that tan^ 0i2 = 0.45 ± 0.05 (3) 

which agrees with the results of the global analysis in table ^ both in the central value and in its 
uncertainty. 

Notice that the only solar model input that enters our approximate determination of solar oscilla¬ 
tion parameters is the solar density around the center of the sun, that controls the 15% correction to 

^SNO measured = (4.9±0.3) 10®/cm^ sec and = (1.74±0.08) 10®/cm^ sec. Each measurement was first 

performed using heavy water (CC and NC events mainly distinguished by their energy spectrum) and later with salt heavy 
water (distinction relies on event isotropy). The measurements performed with these different experimental techniques 
agree. Taking into account the SK measurement of the ES rate, 'l?(r'e) + 0.155'l?(i/,i,.r) = (2.35 ± 0.06) 10®/cm^ sec, the 
SNO measurement of the ES rate, 4>(i^e) + 0.1554>(r',i,T) = (2.36 ± 0.19) 10®/cm^sec and the solar model prediction 
'I>(i^e,(/,T-) = (5.05 ± 0.86) 10®/cm^sec would only marginally improve the measurement of {Pipe —> r'e)) to {P{ve —> 
vfi) = 0.360 ± 0.028. 

^The factor giving the correction to P{ve —» Ve) = shP 612 ranges between 1.1 and 1.2 within the present best-fit 
region at 90% CL. 
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(P(t'e ^ z^e)) in eq. Q. This correction factor is comparable to the Icr uncertainty in {P{ye ^e))'- 
indeed the associated increase of v^) at smaller E^, has not been observed in SNO and SK 

spectra. For the reasons explained above the ‘solar model independent fit’ of m gives now a result 
almost identical to the standard fit, so that we do not show the update of this result. 

Global fits remain still useful for testing if the pieces of data not included in our simplified analysis 
(that have a minor impact in the standard fit) contain statistically significant indications for new 
physics beyond LMA oscillations. At the moment the answer is no. E.g., fig. Eb updates the CPT- 
violating solar fit of m- the best-fit region includes the CPT-conserving limit (diagonal dotted line). 

1.3 Effects of 9 i 3 and low-energy neutrinos 

At low El, matter effects are negligible and the survival probability is given by averaged vacuum 
oscillations. LMA oscillations with 613 = 0 predict the low-energy limit of P{i^e z'e) ia terms of its 
high-energy limit as^ 

P{i'e Vei small Ey) = 1 - 2P{ue z^e, large Ey) + 2P{ue Ue, large Eiy)"^. (4) 

A non zero ^13 as well as new physics beyond neutrino masses allow to avoid this prediction. It is 
therefore interesting to extract these two ideal observables from data. As discussed in the previous 
section P{i'e z'ej large E^) is presently dominantly determined by SNO. The low energy limit of 
P{ve z'e) is presently dominantly determined by Gallium data and can be extracted by a simple 
approximate argument [laiii. Subtracting from the total Gallium rate 

(68.1 ± 3.7) SNU = Rgii = (5) 

its ®B contribution (as directly measured by SNO via GG, = 4.3 ± 1 SNU) and regarding all 
remaining fluxes as low energy ones, suppressed by P(t'e i^ej small Ey), determines it to be 0.57 ± 
0.03. Alternatively, by subtracting also the intermediate-energy CNO and Beryllium fluxes, one gets 
P{i'e —> small El,) = 0.58 ±0.05. We here included only the error on the Gallium rate, which is the 
dominant error. This rough analysis shows that the result only mildly depends on how one deals with 
intermediate energy neutrinos, and on model-dependent details of the intermediate region, thereby 
suggesting the following general result: 

P{iye z^ej large E,^) = 0.31 ± 0.03, P{i^e z/g, small E,^) = 0.58 ± 0.04. ( 6 ) 

As illustrated in fig. El at Ei, = 1 MeV (lOMeV) the difference between the exact LMA prohle of 
P(z^e z^e) and its low energy (high energy) limit is smaller than the la uncertainties in eq. (jHl). High 
energy neutrinos have been dominantly measured by SNO at energies around 10 MeV. 

In order to establish the validity of eq. we compare it with the ‘exact’ result of a global analysis 
of solar data: to perform it we must abandon generality and focus on a specific mechanism that allows 
to avoid the LMA prediction of eq. 0 - We consider the case of a non vanishing ^13 (see also m), 
which gives 


P{l^e- 

z/g, large E,^) 

= sin^ 013 ± cos^ 6*13 sin^ 012, 


(7a) 

P{^e- 

4 Ue, small Ey) 

= = Sin'^ 013 ± COS'^ 013 

1 - i sin^ 2012 

(7b) 


i 


®The transition between the two regimes proceeds at Ei, ~ few MeV, and the high-energy regime is approached 
for El, ^ 20 MeV. At higher energies solar matter effects become comparable to the atmospheric mass splitting, giving 
corrections to solar neutrino rates proportional to 
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Figure 4: The energy-dependent survival prob¬ 
ability predicted by LMA, and how experimental 
data restrict the low-energy and high-energy Urn- 
its of P{Ve l'e,Eu)- 


Figure 5: Best-fit regions at 90, 99% CL (2 dof) 
for Boron and Beryllium fluxes from a global 
solar-model-independent fit of solar data. The 
horizontal ellipse shows solar model predictions. 


In this way the equality in eq. gets replaced by a < inequality. In order to access also the other 
region we analytically continue to imaginary 613 , where cos^ 613 > 1. (This is analogous but less usual 
than allowing sin^ 2023 > 1 in the atmospheric fit). The global fit is performed by keeping fixed Amf 2 
at the central value suggested by KamLAND, because a variation of Am ^2 within the range in tabled 
negligibly affects solar data. This is the key extra input provided by KamLAND; a more complicated 
solar plus KamLAND global fit would give the same result. 

The result of the global solar fit performed in this specific context is shown in fig. |23, and agrees 
with the semi-quantitative general result of eq. ®. We see that the LMA prediction in eq. (® is well 
compatible with data. The constraint on 613 provided by solar data (subdominant with respect to 
the constraint from CHOOZ and atmospheric data) is included in the global analysis summarized in 
table n and fig. Q 

We performed more global fits considering a few other ways of avoiding the LMA prediction of 
eq. (® : in each case the allowed region looks like a potato similar to the one in fig. |2j3. Different fits 
give allowed regions with sizes and shapes that vary roughly as much as potatoes vary. For example we 
tried to linearly distort the i-’(r'e —> t'e, E^) profile predicted by LMA for 613 = 0 as i-’(r'e —> t'e, E^) = 
sin^ 012 + A[PLMA(^'e Vf.,Ey) — sin^0i2]. This is somewhat different than considering a 0i3 0, 

because matter effects depend on 0 i 3 , that therefore does not act as a linear distortion. 

The main point is that the approximate general result of eq. ® fairly summarizes the variety of 
exact results obtained by performing global fits in presence of different mechanism that distort the 
prediction of eq. without introducing new notable features at intermediate energies. Therefore 
eq. ® is a useful semi-quantitative way of summarizing our present knowledge of PiVf. 1 /^). Its 
behavior at intermediate energies 1 MeV ^E^y^lO MeV is basically unknown. 
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1.4 Low energy neutrinos and BOREXINO 


Can future sub-MeV solar neutrino experiments improve on oscillation parameters? This question was 
answered in m, and we do not have much to add; for recent works on the subject see m- We recall 
here the main points. Low energy solar neutrinos do not experience the MSW resonance in the sun: 
their survival probability is therefore given by the averaged vacuum oscillations expression of eq. ()7b|) . 
This means that, in first approximation, sub-MeV experiments have nothing to tell about Am^ 2 ) but 
could give information on 0 i 2 by measuring P{i'e t'ej small Ey). Within the 3 neutrino context 
this same survival probability of eq. (I7bl) is also directly measured by reactor Vg experiments, such as 
KamLAND. KamLAND can achieve a determination of P{ve t'e; small Ey) competitive with the 
solar result in eq. Q. A future reactor experiment with baseline ~ 50km appropriate for observing 
the first oscillation dip could achieve an error (2 3) times lower than the present error in eq. (|H|). 

Such reduced uncertainty is already achieved today, if one trusts the LMA prediction of eq. 0) and 
employs it to infer the low-energy limit of Pi^i'e ^e) from SNO data: 

P(t'e ^ t'ej small Ey) ~ ^ ~ 2 ~ ^ (LMA prediction). (8) 

For simplicity we here assumed 0 i 3 = 0, because its value is not a relevant issue: if 0 i 3 is large enough 
to have a sizable effect, long-baseline experiments will see and measure it so well that only the central 
value of eq. (jHl) (but not its uncertainty) has to be changed. 

The precision of sub-MeV experiments is ultimately limited by the 1% solar model uncertainty 
on the pp flux. Although it presently seems unrealistic to aim at measuring the pp rate with this 
level of accuracy, the above discussion shows that such ultimate precision should be achieved, if sub- 
MeV experiments want to contribute significantly to the determination of 612 ■ Going beyond the 
Sv framework, sub-MeV experiments can make important searches for new physics beyond neutrino 
masses. The reason is that high energy solar neutrinos, detected by SNO, are almost pure 1^2 (tbe 
neutrino mass eigenstate with mass m 2 ), so that low energy experiments are needed to probe new 
physics that dominantly affects z/i. This argument applies e.g. in presence of an extra sterile neutrino 
in wide ranges of its oscillation parameters |23j . 

If LMA is the end of the story the near-future experiment Borexino will not improve the deter¬ 
mination of oscillation parameters and should observe no anomalous day/night nor seasonal variation. 
(KamLAND can also be converted into a solar neutrino experiment). We comment on the impact of 
Borexino (and eventually KamLAND) from the point of view of new physics. As discussed below, 
in presence of generic new physics, existing data poorly constrain the Beryllium rate, so that large 
deviations from the LMA prediction are possible. The rate measured by Borexino can be modified 
by new physics that mostly affects intermediate energy neutrino. Rather than performing a detailed 
analysis of one specific source of new physics, out purpose will be to study this possibility in an approx¬ 
imate but sufficiently general way. To this end, we notice that in Borexino the effects of unspecified 
new physics dominantly manifest as an anomalous flux of Beryllium neutrinos (Borexino also has a 
minor sensitivity to pep neutrinos m)- Therefore, in the next section we study how well present data 
restrict the Beryllium rate. 

1.5 Solar neutrino fluxes 

An eventual deviation of solar neutrino fluxes from solar model LMA predictions could be due to new 
physics in neutrinos beyond LMA, or to (new) physics not included in solar models. For concreteness 
we focus on the second case, and study how well present solar neutrino data determine solar neutrino 
fluxes, updating the results of m- Apparently a detailed global fit seems needed to address this issue. 


but — once again — one can answer to this question in a simple way. Following m the main points 
to consider are the following ones: 

1 . The luminosity constraint allows to precisely predict the pp fluxes. Since no experiment so far 
is sensitive to deviations compatible with the luminosity constraint, we can essentially set pp 
fluxes to their solar model value. 

2. SNO measured the Boron flux. 

3. Only two kind of experiments, Gallium and Chlorine, have measured low-energy neutrino fluxes. 
Therefore, data only constrain two linear combinations of low-energy fluxes. 

4. The Chlorine experiment has a poor sensitivity to low energy neutrinos: after subtracting the 
~ 80% Boron contribution to the Chlorine rate, as directly measured by SNO via CC, the 
residual low-energy contributions to the Chlorine rate is just about 2a above zero. 


Therefore the Chlorine rate carries so little information on low energy fluxes, that our present knowl¬ 
edge on low-energy fluxes is well summarized by a single number: their contribution to the Gallium 
rate. Restarting from eq. we now subtract from the total Gallium rate Rca its ®B contribu¬ 
tion and its pp,pep contributions (as predicted by solar models and LMA oscillations, see eq. ©: 
Rpp^pep = 41.3 ± 1.5), obtaining 


(22.5 ± 4) SNU = ilgfJo + R?i. = 


(9) 


We have taken into account that LMA oscillations suppress both rates by about 0.55 ±0.02 — a value 
negligibly different from the low-energy limit of Pii'e ^e) of eq. (j8|). 

In order to show the accuracy of our simplified analysis, we now compare its results with the 
‘exact’ results of a solar-model-independent global analyses of solar and KamLAND data, performed 
as described in m- To perform such comparison we need to consider a well defined, simple and 
relevant sub-case: we assume a non standard Beryllium flux but a standard CNO flux, $cNO ~ 
0.6 • 10 ^/cm^s. As previously discussed, this choice is motivated by the fact that Beryllium neutrinos 
are more important, that Borexino should study them, and that according to solar models the CNO 
contribution to any measured solar neutrino rates is smaller than its experimental error. In this way 
eq. 0 reduces to $730 = (4.9 ± 1.1) • 10^/cm^s. This can be directly compared with the result of the 
global analysis, shown in fig. El that can be approximatively summarized as 


4>8b = (5.0 ± 0.2) • lOVcm^s, (10a) 

4*736 = (4.1 ± 1.1) • 10 ®/cm^s. ( 10 b) 


The global fit also includes the Chlorine rate, which has a central value about 2 cj lower than the 
LMA prediction, and consequently somewhat reduces 4 > 73 e. From these numbers, we see that the 
determinations of the Boron and Beryllium fluxes that follow from our simplified analysis are quite 
adequate. Borexino will significantly improve over the present determination of 4>73g. 


The simplified analysis also determines our present knowledge of the CNO flux. As previously 
discussed solar neutrino experiments basically measured the linear combination of CNO and Beryllium 
fluxes of eq. ©, that directly implies an upper bound on the CNO flux: 

4*cno < 6 • 10®/cm^s at 3cr (1 dof), (11) 

which is one order of magnitude above solar model predictions. Eq. CU) updates the value first 
obtained in m (section 4). This constraint was also re-derived in |2H] , emphasizing that the result of 
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eq. © or of eq. m proves that the CNO cycle does not give the dominant contribution to the total 
solar luminosity Lq. Indeed by converting the neutrino flux <hcNO into the corresponding energy flux 
eq. CU) reads Lcno ^0.1 Lq. 

2 Oscillations with atmospheric frequency 

Present data do not precisely determine |Am 23 | nor 623 (see table lUor fig. ©, give an upper bound on 
013 and do not determine the sign of (i.e. if neutrinos have ‘normal’ or ‘inverted’ mass hierarchy). 

Several experimental programs using long-baseline and reactor neutrinos plan to to confirm the SK 
evidence and to improve on Am^^, 023 and 0 i 3 , possibly measuring a non-zero value of the latter angle. 
From the point of view of the 3 neutrino framework these experimental programs seem adequately 
complete m, thanks in particular to JHF EH, possibly complemented with a new reactor neutrino 
experiment m Here, we discuss the present determination of atmospheric parameters and other 
implications of the existing data. 

The global fit. The atmospheric fit includes the final SKI data [2], the final Macro data m 
and the latest K2K m data. As in other similar analyses, the SK and K2K fits are extracted from 
the latest SK and K2K papers by ‘graphical reduction’ (i.e. using a scale and a pencil) because this 
procedure guarantees a more accurate result than an independent reanalysis. We employ the ‘standard’ 
analysis of final SKI data. The SK collaboration also performs an alternative analysis, by selecting 
the data which have the best resolution in L/E^,, obtaining similar central values and somewhat 
different uncertainties. (The CPT-violating fit of fig. Of is instead based on an independent reanalysis 
of atmospheric data, so its CPT-conserving limit can slightly differ from the CPT-conserving fit of 

fig-©- 

The parameter 023. The value of this parameter can be obtained from a simple physics argument. 
It is dominantly determined by SK multi-GeV /r-like events as sin^ 2023 — 2(1 — N-^/N^) = 1.02 ±0.08, 
where ~ 400 {N-^ ~ 200) are the number of down-ward (up-ward) going ^ events, that experience 
roughly no oscillation (averaged oscillations). A detailed analysis is needed to include the rest of 
the SK, K2K, Macro data, which however only mildly improve on the determination of sin^ 2023. 
An improved measurement of sin^ 2023 from the up/down ratio of atmospheric neutrinos could be 
performed at a future Mton-scale atmospheric detector.^ Future long-baseline experiments with very 
intense conventional neutrino beams m will also lead to progress. 

The parameter lAm^gl. The present situation concerning |Am| 3 | is quite different: SK cannot 
precisely measure it and cannot see a clear oscillation dip, and a detailed analysis is necessary to extract 
its central value and error. Long-baseline experiments should significantly reduce the uncertainty 
on |A 7 T 7 -| 3 | by identifying the energy at which neutrinos experience the first oscillation dip. This 
measurement has been already performed by K2K, but with poor statistics: K2K achieves m a 
determination of |Am 23 | with central value and error close to the one of SK |2]. 

The parameter 0 i3. The CHOOZ constraint on 0i3 is strongly correlated with the determination 
of |Am| 3 |. As discussed in section ©31 solar data have a subdominant impact on the determination of 
013, comparable to the effect of changing the kind of analysis of SK atmospheric data (we used the 

'^One should take into account the few % effects of solar oscillations and of @13 (see m for a recent discussion) which 
will be respectively better measured by solar and reactor/long-baseline experiments. 
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Figure 6 : 99% CL expected ranges for the parameters rucosmo = fni + m 2 + m 3 probed by cosmology 
(fig.\B(i) andmy^ = probed by jd-decay (fig.\^) as function of the lightest neutrino mass. The 

darker lines show how the ranges would shrink if the present best-fit values of oscillation parameters 
were confirmed with negligible error. 

‘standard’ SK analysis). The statistically insignificant hint for a 6*13 > 0 in fig. □ is mainly due to a 
small deficit of events in CHOOZ data at lowest energies. 

Other effects? Data show no significant hint for new effects beyond three massive neutrinos. For 
example fig. shows a global fit performed without assuming that neutrinos and anti-neutrinos 
have the same atmospheric mass splitting and mixing angle. We see that the best-fit lies close to 
the CPT-conserving limit, and that the atmospheric mass splitting in anti-neutrinos is poorly deter¬ 
mined. Nevertheless, this is enough to strongly disfavor a CPT-violating interpretation of the LSND 
anomaly m- Near-future long-baseline experiments will probably study only z/ rather than n. 


3 Non-oscillation experiments 

In this section we discuss non-oscillation experiments and consider the 3 non-oscillation parameters 
mentioned in the introduction. Making reference to experimental sensitivities, the 3 probes should 
be ordered as follows: cosmology, 0u2f3 and finally /3 decay. Ordering them according to reliability 
would presumably result into the reverse list: cosmological results are based on untested assumptions, 
and Gi'2(3 suffers from severe uncertainties in the nuclear matrix elements. Even more, there is an 
interesting claim that the Qv2f3 transition has been detected m (see section ESI for some remarks), 
there is a persisting anomaly in TROITSK /3 decay, and even in cosmology, there is one (weak) claim 
for a positive effect. None of these hints can be considered as a discovery of neutrino masses. Several 
existing or planned experiments will lead to progress in a few years. 

In this section, we assume three massive Majorana neutrinos and study the ranges of neutrino 
mass signals expected on the basis of oscillation data, updating and extending the results of m- 
Our inferences are summarized in table Q and obtained by marginalizing the full joint probability 
distribution for the oscillation parameters, using the latest results discussed in the previous sections. 


11 























3.1 Cosmology 


Cosmological observations are mostly sensitive to the sum of neutrino masses: rricosmo = mi+m 2 +m 3 , 
that according to standard cosmology controls the energy fraction in non relativistic neutrinos as 
Q-yh? = mcosmo/93.5eV where as usual h parameterizes the present value of the Hubble constant as 
h = lltoday/(100km/sMpc). Cosmology does not distinguish Majorana from Dirac neutrino masses. 

In order to convert CMB and LSS data into a constraint on neutrino masses one needs to assume a 
cosmological model. The “WMAP constraint” m assumes that the observed structures are generated 
by Gaussian adiabatic primordial scalar fluctuations with flat spectral index n evolved in presence of 
the known SM particles, of cold dark matter and of a cosmological constant. This standard model of 
cosmology seems consistent with all observations. Furthermore, it is assumed that observed luminous 
matter tracks the dark matter density up to a scale-independent bias factor. The data-set includes 
Lyman-a data about structures on scales so small that perturbations are no longer a minor correction 
to a uniform background. These data are sensitive to neutrino masses (and thereby somewhat affect 
the global cosmological fit) but might be affected by non linear evolution effects, which are not fully 
understood. In summary, cosmology presently gives the dominant constraint, which however rests on 
untested assumptions and on risky systematics. 

In the future the sensitivity to neutrino oscillations will improve thanks to better CMB data and 
to new LSS measurements less plagued by potential systematic effects. If cosmology were simple (e.g. 
a spectral index n = 1 , no tensor fluctuations,...) then it seems possible to detect even neutrino 
masses as small as allowed by oscillation data j32]- The expected ranges of mcosmo are reported in the 
lowest row of tabled in the limiting case where the lightest neutrino is massless, and in fig. EK in the 
general case. 

3.2 Direct search via (3 decay 

/3-decay experiments dominantly probe the quantity rriy^ = (m • . If neutrinos are quasi¬ 

degenerate, rriy^ is their common mass. The MAINZ and TROITSK experiments obtained comparable 
limits: = —1.2±2.2±2.1 eV^ jS3] and = —2.3±2.5±2.0 eV^ j31] respectively. The 95% bounds 

quoted by the experimental collaborations agree with the values obtained in Gaussian approximation. 
Thus, we combine the two measurements by summing errors in quadrature and get 

ml^ = -1.7 ± 2.2 eV^ i.e. rriy^ < 2.0 eV at 99 % C.L. (12) 

In order to study what oscillation data imply on the value of rriy^ we write it in terms of the neutrino 
masses and of the mixing angles 9ij as 


(13) 


In the case of normal mass hierarchy, mi <C m 2 m 3 oscillation data imply the 99% CL range 
niy^ = (4.6-^ 10) meV. In the case of inverted hierarchy, m 3 ^ mi < m 2 one gets niy^ = (42-^57) meV. 
The last number if a factor 5 below the planned sensitivity of Katrin [SSI. 

It is immediate to obtain the ranges corresponding to the generic case of a non vanishing lightest 
neutrino mass mughtest (where mughtest = 'mi in the case of normal hierarchy and mughtest = m 3 in the 
case of inverted hierarchy). As clear from the definition m^^ = (m • m^)ee or from the more explicit 
expression in eq. m one just needs to add w-fightest resulting bands at 99% C.L. are 

plotted in fig. IBJa. 
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Nucleus and 
experiment 

observed 
events, n 

background 
events, b 

expected 

signal 

99% C.L. bound 
on \mee\/h 

^®Ge HM [SSI 

21 

20.4 ± 1.6 

76 mee/eVp/h^ 

0.44 eV 

7®Ge IGEX jSnillll 

9.6 

17.2 ±2 

23.5 mee/eVp//i^ 

0.55 eV 

^®°Te CUORICINO dSlEZl 

24 

35.2 ±4 

21.5 mee/eVp/h^ 

0.62 eV 


Table 2: Numbers of observed events, expected background and predicted signal in the most sensi¬ 
tive Onip experiments. The last column shows the 99% C.L. eonstraint on \mee\/h, where /i ~ 1 
parameterizes the uneertain OvljS nuelear matrix element and depends on the nucleus studied. 

3.3 Neutrino-less double-beta decay 

Updating the results of 1301, in the 3 neutrino framework we discuss the connection with oscillations, 
the bound from 0i^2/3 on neutrino masses and the possible hint of a signal. 

IGEX, CUORICINO and nuclear uncertainties. First, we recall that recently two experiments 
produced new relevant data. The first is IGEX (SE], now terminated, which used 86 % enriched ^®Ge 
(with exposure 7 • 10^^ nuclei-yr). The second is CuORiCiNO EH, now running, which uses 34% 
natural ^^°Te (with exposure 1.4-10^'^ nuclei-yr, and the plan to collect 10 times more data in about 6 
years). Both experiments reached a good level of background (about 0.18/keV-kg-yr), good detection 
efficiencies of 70% and 84%, and respectable energy resolutions, the full width at half maximum 
being 4 keV and 7 keV respectively. For comparison, the Heidelberg-Moscow (HM) data-set with 
pulse-shape-discrimination m has exposure 2.5 • 10^® nuclei-yr and background 0.06/(keV - kg - yr). 

As in 1301 we quote the limits on |mee| adopting the 0z^2/3 nuclear matrix elements A^o computed 
in 1301- To use a different calculation with matrix element A4, just rescale by the factor h = MojAd, 
which depends on the nucleus studied.® We always explicit the factors h when quoting an experimental 
result on t)v2(5. We prefer to show such uncertainty explicitly rather than attempting to evaluate the 
theoretical error on matrix elements.® Different published computations find h in the following ranges 

0.3 gni < h(™Ge) < 2.4 dH, 0.4 gOl < /i(^®°Te) < 2.7 (14) 

When comparing two results on \mee\ obtained with different nuclei one needs to consider the ratios 
of h, e.g., /i(^®®Te)//i(^®Ge). This quantity is also uncertain, spanning the following range: 

0.3 HE] < /i(i®0Te)//i(^®Ge) < 1.7 gg. (15) 

Now we convert Qu2f3 data into a constraint on \mee\ following the same simple procedure employed 
by the HM collaboration I3H1: the Poisson likelihood of having s signal events is C{s) oc e ^{b s)"', 

where n and b are the numbers of observed and expected background events in the Scte window around 
the Q value of the t)v2(5, and ue is the energy resolution of the apparatus. We therefore evaluate n 
and b in the 10 keV region around Q = 2039 keV for IGEX |36[ I44| and in the 18 keV region around 

Q = 2529 keV for CuORiCiNO gEllSZl- Results are collected in tabled For both experiments the 

observed number of events is slightly below the expected background. In this paper we systematically 
employ the Gaussian approximation: the 99% constraint (2.58 (t) on \mee\ is computed by defining 

= — 21 n/l, marginalizing it with respect to the uncertainty in the background, and finally quoting 
the value of |mee| such that x^(mee) — = 0) = 2.58^. Different approaches to statistical 

inference give slightly different constraints. 

®The symbol h has different meaning for cosmology and for Oi/2/3, as should be clear from the context. 

®It would be useful if future calculations of matrix elements could provide error estimates. 
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Since HM and IGEX are both based on ™Ge, it is possible to obtain a more stringent combined 
bound, \mee\ < 0.38 h eV7 

Inference on |mee| from oscillations. Using the latest oscillation data, we study the expected 
range of |mee|. We follow the procedure described in m- Assuming three CPT-invariant massive 
Majorana neutrinos we get: 

a) In the case of normal hierarchy (i.e. mi -C m 2 ^ m 3 , or > 0) the ee element 

of the neutrino mass matrix probed by Oz/2/3 decay experiments can be written as \mee\ = 
|g 2 iOj^sun _|_ g 2 i/ 3 j.^atm| ^ g^j-g unknown Majoi’ana phases and the ‘solar’ and ‘atmo¬ 

spheric’ contributions can be predicted from oscillation data. The solar contribution to mee 
is sin^ 012 cos^ 013 = (2.78 ± 0.22) meV. The bound on 0 i 3 implies rUgg™ = 

sin^ 0 i 3 < 1.7meV at 99% CL. By combining these two contributions we get the range 

reported in tabled The precise prediction is shown in fig.C^. Unlike in our previous results m, 
a non zero contribution is now guaranteed at high confidence level, because data now tell that 
the ‘solar’ contribution is larger than the ‘atmospheric’ contribution, so that a cancellation is 
not possible. 

b) In the case of inverted hierarchy (i.e. m 3 mi ~ m 2 , or Am^g < 0), from the prediction 

\mee\ ~ |Am|g|^/^ X I cos^ 012 -|-sin^ 012 1 cos^ 013 we get the range reported in table Q (see 

also m)- The precise prediction is shown in fig. [33. The main uncertainty on |meg| is due to 

the Majorana phase a, rather than to the oscillation parameters. If the present central values 
were conhrmed with infinite precision we would still have the loose range |mee| = (19 50) meV 

at any C.L. 

The I mee I ranges for normal and inverted hierarchy do not overlap. Values of |mee| outside these 
ranges are possible if the lightest neutrino mass is not negligible, as shown in fig. |3i. In the inverted 
hierarchy case (Am|g < 0) a non zero lightest neutrino mass can only make |mge| larger than in case 
b). The darker regions in fig. |3i show how the predicted range of |mge| would shrink if the present 
best-fit values of oscillation parameters were confirmed with infinite precision. The two funnels present 
for Am 23 > 0 have an infinitesimal width because we assume 0 i 3 = 0 and would have a finite width 
if 013 / 0. 


c) The case of quasi-degenerate neutrinos with mass m^^ corresponds to the upper region of 
fig. 131 . |meg| and rngosmo are related to m^^ as 

rucosmo = 3mi/g and 0.24 mi/^ < |mge| < m^^ at 99% C.L. (16) 

The lower bound on |mge| holds thanks to the fact that solar data exclude a maximal solar 
mixing and that CHOOZ requires a small 0i3. 


Upper bound on neutrino masses from 0u2f3. The above result means that by combining 
oscillation data with the Gi'2[3 upper bound on |mee| implies an upper bound on the parameter m^,^ 
probed by (5 decay experiments m- In view of present values, such m^,^ corresponds to the common 
mass of quasi-degenerate neutrinos. This bound is shown in fig.|7[:, and depends on the 0u2f3 nuclear 
matrix elements, parameterized by the uncertain factors /i ~ 1. Therefore the results of CUORICINO, 
obtained with a nucleus different than HM and IGEX, add confidence in the result. For h = 1 the 

cannot combine HM/IGEX with CuORiCiNO in a reliable manner, due to the uncertainty on the relative nuclear 
matrix element discussed above. Using the matrix elements of |.‘-{h| for ^®Ge and ^'^°Te, would give \mee\ < 0.34 h eV at 
99% G.L. with h = 1. 
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Figure 7: Predictions for |mee| assuming a hierarchical (fig. w and inverted (fig. w neutrino spec¬ 
trum. In fig. we update the upper bound on the mass of quasi-degenerate neutrinos implied by 0u2f3 
searches. The factor /i ~ 1 parameterizes the uncertainty in the nuclear matrix element (see sect. 
2.1). In fig. we plot the 99% CL range for mee as function of the lightest neutrino mass, thereby 
covering all spectra. The darker regions show how the m^e range would shrink if the present best-fit 
values of oscillation parameters were confirmed with negligible error. 
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Figure 8; Fig. 03: the latest HM data m f’Tl.Tkg • yr) used to claim a 4.2cj evidence for 0u2(3. 
Fig. 0); the statistical significance of the t)v2f5 signal, as function of the assumed flat component of 
the background. 

combined constraint is my^/h < 1.0eV at 99% C.L. This constraint is stronger than the /3-decay 
constraint (but holds under the additional assumption that neutrinos are Majorana particles) and is 
weaker than the cosmological constraint (but needs no assumptions about cosmology). 

In all above discussions we employed the bound on |mee| from HM data, as published by the HM 
collaboration m- 

Remarks on the hint for Oz/2/3. While the HM collaboration used their data to set a bound on 
|mee|, some members of the HM collaboration reinterpreted the data as an evidence for 0u2f3 m- 
According to this claim, the latest HM data m plotted in fig. |HK contain a 4.2(T evidence for a Qv2/3 
peak (indicated by the arrow). In these latest results the peak is more visible than in latest published 
HM data 131, partly thanks to higher statistics (increased from 53.9 to 71.7 kg yr) and partly thanks 
to an ‘improved analysis ’ m)- 

This claim is controversial, mainly because one needs to fully understand the background before 
being confident that a signal has been seen. In order to allow a better focus on this key issue, we 
present fig. ISb) which should be uncontroversial. It shows the statistical significance® of the 0i^2/3 
signal as a function of the true background level b, assumed to be quasi-flat close to the Q-value of 
0i'2f3, Q ~ 2039keV. The crucial issue is: how large is b? The HM collaboration earlier claimed |38| 
b = (13.6 ± 0.7)events/(71.7kg yr • keV). In such a case the statistical significance of the signal would 
be less than Icr, see fig.|Eb- This can be considered as the upper bound on b computed assuming that 
all events in a wider range around Q come from a quasi-flat background. 

A statistically significant hint for 0z^2/3 is obtained if one can show that b is lower. The continuous 
line in fig. 0 shows a fit of HM data using a tentative model of the background assumed to have 
a quasi-flat component b (mainly due to ‘natural’ and ‘cosmogenic’ radioactivity) plus some peaks due 
to faint 7 lines of which is a radioactive impurity present in the apparatus (from the 

chain). Their positions and intensities can be estimated from tables of nuclear decays; however they 
are modified by 0(1) factors by detector-related effects which depend on the unknown localization of 
(see m and ref.s therein). The fit in fig. 0 is performed by allowing the intensity of each line 
to freely vary. In this way part of the background is interpreted as ^^^Bi peaks, thereby reducing the 

^Defined in Gaussian approximation as (Xno signal “Xbest where = —21n£ and C is the likelihood computed 

combining statistical Poissonian uncertainties with other systematic uncertainties. 


16 















































quasi-flat component b. Proceeding in this way, we find that the statistical significance of the 0u2f3 
signal is about 2.7a. This is the best we can do with the data at our disposal. This kind of analysis 
was proposed in m and has been adopted in m- 

However, some details in its implementation prevent this analysis from fully reaching its goal, 
which is determining b from regions with no peaks. 1) The latest data have been published only 
below 2060 keV. (Above 2060 keV in fig. Et- we plotted HM data artificially rescaled to account for 
the larger statistics). Below 2060 keV there are little energy ranges with no peaks. Including the old 
data above 2060 keV in the fit would reduce the significance of the signal down to about 2.2cr. 2) HM 
data contain hints of extra unidentified spurious peaks at specihc energies (at 2030 keV and above 
2060 keV). Fitting data assuming that these extra peaks can be present at arbitrary energies with 
arbitrary intensities reduces b and enhances the statistical significance of the signal. 

Various future experiments plan to test the claim of m- The most direct test requires using the 
same technique (germanium detectors) but in the worst case, when the observed hint of a signal is 
due to some irreducible (hypothetical) background, a safe test requires a different type of detector. 
As previously discussed CuORiCiNO employs a different nucleus: its capabilities relative to HM and 
IGEX depend on the uncertain relative nuclear matrix element. For instance, with the lowest value in 
eq. da CuORiCiNO is already testing the claim of m-, for intermediate values new data of CuORiCiNO 
will significantly test the claim; for the highest value CuORiCiNO will be not sufficient. 

In our view, a discussion of what should be considered as a convincing evidence for Oz/2/3, is anyway 
useful or necessary, because any experiment (past and future) needs to confront with this issue. 


4 Conclusions 


Assuming oscillations of three active massive neutrinos we updated the determination of the oscillations 
parameters, at the light of latest experimental data. Results are shown in tableland in fig. ^ We 
notice that the parameters Am^ 2 ) ^12 > ^23 are now dominantly determined by simple and robust sub¬ 
sets of data, such that simple arguments give the same final result as global analyses. Pieces of data 
that play a sub-dominant role in parameter determination allow to test our assumptions: e.g. allowing 
neutrinos and anti-neutrinos to have different masses and mixings gives the CPT-violating fit of fig. |31 
Present data do not contain evidence for extra effects. 

We also tried to study in a simple and general way related topics, such as the determination of 
the limiting survival probability of solar neutrinos at small and large energies, the present knowledge 
of solar neutrino fluxes, etc. Approximated general results have been compared with ‘exact’ results of 
global fits performed in specific cases. 

Finally, updating the results of m, we studied how oscillation data allow to infer the combination 
of neutrino masses probed by cosmology, /9-decay and 0i^2/9-decay experiments, and discussed the 
present experimental situation. 
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